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ABSTRACT 

When a controller is implemented in a digital computer, with A/D and D/A conversion, 
the numerical errors of the computation can drastically affect the performance of the control 
system. There exists realizations of a given controller transfer function yielding arbitrarily large 
effects from computational errors. Since, in general, there is no upper bound, it is important to 
have a systematic way of reducing these effects. Optimum controller designs are developed 
which take account of the digital round-off errors in the controller implementation and in the 
A/D and D/A converters. These results provide a natural extension to the LQG theory' since they 
reduce to the standard LQG controller when infinite precision computation is used. But for finite 
precision the separation principle does not hold. 


91 



I. INTRODUCTION 


LQG controllers are normally designed under the assumption that computer implemention 
will be perfect (this is the infinite wordlength assumption for state variable computation). 
However, real control systems are subject to the effects of finite wordlength computation. These 
round-off errors should not be ignored in the design of the controller. The influence of these 
errors on the control system and the optimum controller design considering their effects are the 
subjects of this paper. 

We consider the problems that arise with fixed-point arithmetic and the finite word length 
of digital computers. This paper was motivated by the work of Kladiman and Williamson 
[1989]. Mullis and Roberts [1976] and Hwang [1977] in the field of signed processing first 
revealed the fact that the influence of round-off errors on digital filter performance depends on 
the realization chosen for the filter implementation. To minimize round-off errors these papers 
suggest a special coordinate transformation T prior to filter (or controller) synthesis. 

This is in stark contrast to frequency domain approaches to control, which regard as 
irrelevant (and hence is completely ignored) the state space realization of the controller transfer 
function. 

The idea of applying a coordinate transformation prior to controller synthesis has been 
applied to Kalman filter and LQG controller design problems, Williamson [1985], Kladiman and 
Williamson [1989]. One may select the wordlength of the computer to insure that the resulting 
degradation in the performance from round-off error is less than a certain percentage of the ideal 
behavior of the standard Kalman filter or LQG controller without round-off error. This approach 
was adapted by Sripad [1981] in the design of Kalman filters, and later by Moroney, et. al [1983] 
for LQG controller design. In these papers the standard Riccati equations are solved, followed 
by a coordinate transformation to reduce the effects of round-off errors. We shall call these 
controllers LQGx to indicate a standard LQG controller followed by an "optimal" coordinate 
transformation T. This transformation depends on the control gains, hence, we put the word 
optimal above in quotes, because the standard LQG gain is not the optimal gain for the round-off 

92 


error problem. The optimum solution is to design the controller which directly takes into 
account the round-off errors associated with a finite word length implementation, rather than 
merely performing a coordinate transformation T on the LQG controller after it is designed. The 
optimal state estimation problem was solved by Williamson [1985]. This leads to a modified 
Kalman filter. The problem of optimum LQG controller design in the presence of round-off 
error was studied by Kadiman and Williamson [1989], This paper worked with upper bounds 
and numerical results showed improvement over earlier work, but their algorithm does not 
provide the necessary conditions for an optimal solution. This paper provides the necessary 
conditions and a controller design algorithm for the solution of this problem. We shall call this 
controller LQGpw • 

With a fixed point implementation, the states of the LQGpw controller are properly scaled 
to reduce the possibility of overflow. There are many scaling criteria available. The method we 
shall use is the variance oriented procedure, ^-norm scaling [Hwang 1977]. We assume round- 
off errors are additive. This tends to be supported by the literature on state quantization, whereas 
quantization of coefficients leads to multiplicative errors [Williamson 1985], 

The organization of the paper is as follows. In Section 2, the problem of LQG controller 
design in the presence of round-off errors is formulated. The importance of the coordinates of 
the controller will be discussed in Section 3. Section 3 summarizes the needed results from 
[Kadiman and Williamson 1989], and our new results on upperbounds of finite wordlength 
effects. It is shown that the portion of the LQG cost contributed by these errors will range from 
arbitrarily large to an achievable lower bound with the variation of the realization of the 
controller (variation of the choice of coordinates). The coordinate achieving the lower bound is 
described. In Section 4, the optimization problem is discussed in terms of chosing both the 
controller parameter matrices and the realization coordinate simultaneously. The necessary 
conditions are derived for the optimization problem. An algorithm is then presented for the 
designs of the optimal LQGpw controller. The standard LQG and the LQGpw controller are 
compared in Section 5. Some conclusions appear in Section 6. 
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II. Round-Off Error and LQG Controller Design Problem 


In this section, we formulate the LQG controller design problem when round-off errors are 
present. The formulation procedure follows the original ideas of Mullis [1976], Hwang [1977] 
and the ideas of Williamson [1985], Kadiman and Williamson [1989]. Let us assume, for the 
study of round-off error, the discrete controller is designed from a discrete model of the plant to 
be controlled. We then introduce a model for finite wordlength effects into the discrete design 
problem. 

Considering the following discrete-time model of a time-invariant plant: 

Xp(k+1) = ApXp(k) + BpU(k) + DpW p (k) 

<Zp(k) = MpXp(k) + v p (k) (!) 

y p (k) = CpXp(k) 

where x p is the state n p — vector, u, y p and Zp are the control n u — vector, output ny— vector, 
measurement n z — vector, v p and w p are assumed to be mutually independent, zero mean, discrete 
white Gaussian noises with covariance matrices V p and W p , respectively. 

The controller that one might desire to implement is described by following equations: 

jxc(k+l) = A c x c (k) + B c Zp(k) 

{ u(k) = C c x c (k) + D c Zp(k) (2) 

where x c is the controller state n c — vector, u and Zp are the control and measurement vectors 
described in the plant model. In a finite wordlength digital computer, the controller state x c and 
measurement variable will be quantized at each time of computation. Considering the 
quantization process, computation (1) and (2) cannot be accomplished. Instead the computation 
is described by 
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x p (k + 1) = A p Xp(k) + B p Q[u(k)] + DpWp(k) 
Zp(k) = MpXp(k) + v p (k) 
y p (k) = CpXp(k) 


(3a) 


Jx c (k + 1) = A c Q[x c (k)] + B c Q[zp(k)] 

\u(k) = C c Q[x c (k)] + D c Q[zp(k)] (3b) 

where Q[*] stands for the quantization process. Assuming an additive property of the round-off 
error, we can model the quantization process by: 


Q[u(k)] = u(k) + e u (k) 

D/A 

(4a) 

Q[x c (k)] = x c (k) + e x (k) 

control computer 

(4b) 

QtZp(k)] = Zp(k) + e z (k) 

A/D 

(4c) 


where is the round-off error resulting from D/A conversion, e x (k) is the error resulting from 
quantization and e z (k) is the error resulting from A/D conversion. We do not claim that this 
assumption is always justified, but we invoke this common assumption in this paper, since one 
cannot optimize with respect to coefficient errors directly. One can only evaluate designs with 
respect to coefficient errors. There are many such evaluations in filter theory, and we shall add 
our own numerical evaluation in this paper. All such evidence points to a conclusion that 
controller structures that are good with respect to state quantization tend to also be good with 
respect to coefficient quantization. 

It was shown [Sripad 1977] that, under sufficient excitation conditions, the round-off error 
e x (k) can be modeled as a zero mean, white noise independent of w p (k) and v p (k), with 
covariance matrix E x , 


Ex = ql, 



(5a) 


where (3 is the wordlength of the control computer. Similarly, we assume the D/A conversion 
error e u (k) and the A/D conversion error e z (k) to be zero mean, mutually independent white 
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noise and also independent of w p (k), v p (k) and e x (k) with covariance matrices E u and E z , 


E u QuE 

4j_ 2 - 2 P« 
Qu 12 

(5b) 

Ez ~ Qz 

q z 4J- 2 " 2P * 
Hz 12 

(5c) 


where fl u and p z are the wordlengths of D/A and A/D converters. Substitute (4) into (3) to 
obtain a closed-loop system model including finite wordlength effects, 

x p (k + 1) = A p x p (k) + B p u(k) + D p w p (k) + B p e u (k) 
i Zp(k) = M p x p (k) + v p (k) 
y p (k) = CpXp(k) 

x c (k + 1) = A c x c (k) + B c z p (k) + A c e x (k) + B c e z (k) 
u(k) = C c x c (k) + D c Zp(k) + C c e x (k) + D c e z (k) 

We seek the controller to minimize the following cost function 

J = lim E { y p (k)Q p y p (k) + u*(k)Ru(k)} (7) 


(6a) 


(6b) 


where u and y p are again control and output vectors, and Qp and R are the weighting matrices. 

After combining (6a) and (6b), and using the following notation for the vectors and 
matrices: 


x(k) = 


D = 


M = 


x p (k) 

x c (k) 


; y(k) = 


y P 00 

u(k) 


; a = 


Ap 0 
0 0 


, B = 


Bp 

0 


, c= 


Cp 

0 


ii 

O 

i 

D c C c ' 
B c A c 

II 

o 

►— 4 

0 o' 

1 0 

;ii = 

i 

0 

II 

o’ 

I 

J l 
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r 

1 1 

1 




i i 

L _ 


D„ 


M p 0 


Q = 


Q P 

o 


o 

R 


o' 

0 



the closed-loop system is compactly described by 



x(k+ 1) = [A + BGM]x(k) + Dw p (k) + BGIiv p (k) + BGI 2 e x (k) + BGI]e z (k) + Bl^Ck) 

y(k) = [C + IoGM]x(k) + I 0 GIi v p (k) + I 0 GI 2 e x (k) + I 0 GIie z (k) 

and the cost function (7) may be written 

J= lim£{y*(k)Qy(k)} . (10) 

k— 

Now, substitute (9) into (10), since e u (k), e x (k), e z (k), w p (k), and v p (k) are mutually 
independent, 

J = tr{X[C + I 0 GM]‘Q[C + I 0 GM] } + tr{ VpdoGIO'QdoGl!)) 

+ tr{E x (IoGI 2 ) QI 0 GI 2 ))+tr{E z (l o GI 1 ) Q(I 0 GIj ) } (11^) 

where X is the state covariance satisfying: 

X = [A + BGM]X[A + BGM]* + DW p D* + (BGIi)V p (BGIi) 

+ (BGI 2 )E X (BGI 2 )* +(BGI 1 )E z (BGI 1 )*+BI 1 E u (BI 1 )* (lib) 

We can decomposite J in eqn. (11a) into two terms: 

J = J WV + J £ ( i2a ) 


where 

^{XjtC + IoGMf Q[C + 1 0 GM]} +tr{(V p +E z )(I 0 GIi)*Q(I 0 GI 1 )) (12b) 

Xj = [A + BGMJXj [A + BGM]* + DW p D* + (BGIi)(V p +E z )(BGIi)‘ +BI 1 E U (BI 1 ) (12c) 

and 

J e 4 tr(X e [C + IoGM]*Q[C + IqGM]} + tr{E x (IoGI 2 )*Q(IoGI 2 )} (12d) 

X e = [A + BGM]X e [A + BGM]* + (BGI 2 )E X (BGI 2 )* (12e) 

where X = X] +X e . Jwv is the portion of the performance index contributed by disturbances 
eu(k), e z (k), w p (k) and v p (k). J e is the portion contributed solely by round-off error e x (k). 
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To prevent the overflow in controller state variable computation, we must properly scale 
the state variables. We use the ^-norm scaling procedure which is written as: 

[Xi(2,2)]ii = s i = 1, ••• , n c ( 13 ) 

where X, (2, 2) is the (2.2) subblock matrix of X : matrix (the controller subblock), and [-] u 
stands for the ith diagonal element of the matrix. Equation (13) requires that the controller state 
variables have variance equal to s when the closed-loop system is excited only by outside 
disturbance and measurement noise. We call (13) the scaling constraint. 

Therefore, the optimization problem is 

min J = min (J^ + J e ) , ( 14) 

G G v J 

subject to (12-13). 
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ELI. Contribution of Round-Off Error to the LQG Performance Index 


In this section, we discuss the J e term in (12a) and defined by eqn. (12d) which is the 
portion of the LQG cost function contributed by round-off errors. This portion of the cost 
function is coordinate dependent. It is unbounded from above, (that is, it can be arbitrarily 
large), but it has an achievable lower bound, which can be achieved in an optimal coordinate. 
The lower bound result was obtained by [Moroney eL al. 1983] and [Kadiman and Williamson 
1989]. The construction of this optimal coordinate is discussed in this section, where we assume 
G is some given matrix (we shall optimize G later). 

We will first present three key lemmas, which form the basis for the results of this section. 
Lemma 1 . [Mullis and Roberts 1976, Hwang 1977]. 

Given any nxn matrix M, there exist a ( non-unique ) unitary matrix U such that (UMU*)jj = s 
for all j, if and only if tr(M) = sn 

□ 


Lemma 2 . [well known] 

For any two positive definite matrices P and Q, let >4 [•] denote the i 1 * 1 eigenvalue of matrix [•]. 
Then, 

a) MQP] > 0 for all i 

b) The L,[QP] are invariant under the transformation P = TPT* and Q = T~*QT~ 1 where T is 
nonsingular . 

□ 


Lemma 3. 


Let a scalar J be defined by 



J = tr{TT*P} 


(15a) 


where the n p x n p nonsingular matrix T is constrained by 

(T” 1 1~*)ij = s for all i (15b) 


and P is a positive definite matrix. Then over the set of all nonsingular matrices T constrained 
by (15b), 

a) J is not bounded from above. 

b) J is bounded from below (J > J) by 


l =— M ^)] 2 

sn p 


where 


p=4 p Vp 


and Vp" is symmetric. 


c) J in (16a) is achievable by the matrix T; 


T = T = U t n,V* 


where U t , V t are unitary, n t diagonal, satisfying 


u t nr 2 u; 


snp'N/P 

tr(^P) 


(16a) 


(16b) 


(17a) 


(17b) 


[Vtn^V*]^ for alii . 


(17c) 

□ 
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Statements b) and c) are minor modifications of the results obtained by [Mullis and 
Roberts 1976] and [Hwang, 1977]. The proof of a) appears in Appendix A. An algorithm for 
solving (17b), (17c) is given in Appendix B. 


The contribution of finite wordlength error in the cost function is described by equations 
(12d) and (12e). This J e term can also be written as: 

J e =tr{Ke(BGI 2 )E x (BGI 2 )*} +tr{E x (IoGI 2 )*Q(IoGl 2 )} (18a) 

K e = [A + BGM]*K e [A + BGM] + [C + I 0 GM]*Q[C + I 0 GM] . (18b) 

Since E x = ql, we then have: 

J e = qtr{(BGI 2 )*K e (BGI 2 ) + (I 0 GI 2 )‘QaoGI 2 )} . (19) 

We can easily check that the (2, 2)th subblock matrix of Kg (the controller subblock Kg (2, 2)) 
satisfies: 

Kg (2, 2) = (BGI 2 )*K e (BGI 2 ) + (I 0 GI 2 )*Q(I 0 GI 2 ) . (20) 

Substituting (20) into (19) reduces (19) to 

J e = qtr[Kg(2, 2)] . 

Hence, the minimization of J e reduces to the problem: 

minJ e , J e = qtr{K e (2, 2)) (21) 

subject to (18b), (13) and (12c). From the singular value decompositions 

X 1 (2,2) = UjL x U x (22a) 

L^U x Ke(2, 2)U*Zjf = U£2*U k (22b) 

then U x , U k are unitary, L x , £ k are diagonal and 

101 



(22c) 


Xk =diag {... ?ti[K e (2, 2)X,(2, 2)] ... ) . 


Suppose we begin our study with the closed-loop coordinate transformation T as: 


T = 


I 0 

0 u;s*Uk 


(23) 


Then, after this coordinate transformation as suggested by Kadiman and Williamson [1989]: 

Xi(2, 2) = (U:z‘ A Uk)- 1 X 1 (2, 2)(U*X^Uk*r* =1 (24) 

Ke(2, 2) = (UxX^Uk)X(2, 2)(UxX^Uk) = X* . (25) 


If we take one more controller coordinate transformation T c , the index J e and its constraint 
equations, (after we substitute (24) and (25) into (13) and (21)), become 

J e = qtr[T c T*X k ] (26a) 

[T c T c s, i — 1, ■ , n c . (26b) 

Since, from Lemma 2, in (22c) is coordinate independent, we may ignore the K<. and X] 
calculations (18b) and (12c) and concentrate on T c in (26). Then, by applying Lemma 3 on 
equation (26), we have following theorem. 

Theorem 1 . The round-off error term J e in the LQG performance index (12d) and (12e), and 
constrained by the scaling constraint eqn. (12c), (13), is controller coordinate dependent. It is 
unbounded from above when the realization coordinate varies arbitrarily. It is 
bounded from below by the following lower bound: 


J 

“e 



(27) 


The lower bound is achieved by the following controller coordinate transformation: 


t=u;z,xuav; 


(28a) 


where U x , U k , U t , V t are unitary matrices, X x , II, are diagonal matrices, subject to the 
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constraints : 


Xi(2, 2) = U;i x U x 

(28b) 

X^U x Ke(2, 2)U x X)f = Uk^Ujc 

(28c) 

_ SI1 C Xfc 

u ' n - u '= s* 

(28d) 

[V,nr 2 vrfc = s. i=l. ••• ,nc ■ 

(28e) 


□ 


To find the optimal coordinate transformation T. in (28a), we must solve (28d), (28e) to 
obtain U t , n t , V t . The equations (28d), (28e) are, however, special cases of (17b), (17c), where 
P is the diagonal matrix X k . An algorithm is given in Appendix B to compute the U t , n t , V t 
needed for (28a). 

The conclusion of this section is that the problem min J e is solved by the coordinate 

T c 

transformation given by (28a). 
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IV . LQG Controller Design in the Presence of Round-Off Errors 


As discussed in Section n, when round-off error is present, the LQG performance index 
can be decomposed into two terms. One term contains the influence of disturbance and 
measurement noise, the other term is contributed by round-off errors. Although the first term is 
not influenced by the coordinate of the controller, the second term is critically dependent on the 
coordinate. An optimal coordinate transformation is given by (28a). With the scaling 
requirement of the controller state variables to prevent overflow, we have a different 
optimization problem now for controller design comparing to the original optimal control design 
problem without round-off errors. In this section, we will discuss the controller design. 

Let us first present a useful result. 

Lemma 4 . Suppose = £ -^[ K(i, i)X(j, j)] where K(i, i) and X(j, j) are the (i, i)th subblock 
of K and (j, j)th subblock ofX respectively. Define 

then: 


a ) ^k J kx(P.q) = 0 when p * i or q * i 


(29a) 


VJ fa n\- 1 y IE ] (i, j)]/th -row [E(i, jM/VccXO, j) 
kJkxtp.q; - 7 L P-- ■ when p = i and q = i (29b) 

-film, i)x(j, j» 


b) Jkx (P> Q) = 0 when p * j or q * j 


(29c) 


X7 T. nU ^ v K(i, i)[E (i, j)]/th-row[E(i, j)]/th-col 

v xJkxtP, qi- — 2, ;= = When p=jandq=j (29d) 


' i = l 


•\/X/[K(i, i)X(j, j)] 


where V k J kx (p, q) and V x J kx (p, q) are the (p, q)th subblock ofV k and V x J kx , E(i, j) is the 
eigenvector matrix of matrix K(i, i)X(j, j) 


J04- 


□ 


The proof of the lemma is given in Appendix A. 

The LQG controller design problem, when finite wordlength effects are taking into 
account, are described by the equations (12-14). This is denoted as the LQGpw controller. 
However, the scaling constraint (13) can be always satisfied by properly choosing the 
coordinates of the controller, so the problem breaks up into two parts: Finding G and finding its 
optimal coordinate transformation T c to satisfy (12), (13) and (14). On the strength of Section 3, 
we can therefore write the optimization problem as 

min J = min (Jwv + J e ) = min[min (J^ + J e )] 

G,T C G,T C G T c 

since is constant in terms of the variation of T c , we have 

min J = min [J^ + min J e ] (30) 

G,T C G T c 

Assume J =minJ e is given by (27), from Theorem 1. Hence, the equivalent LQGpw design 
T c 

problem becomes 


minfJwv + J ] , (30a) 

G ~ t 


subject to (12c) and (18b) where 

3^ = trX, (C + I 0 GM)*Q(C + I 0 GM) + tr(V p + EJOoGIj j'QaoGIj ) (30b) 

J = — (trlk) 2 (30c) 

^ Site 

where L k is defined by (22c), and the transformation T c which yields is given by the 

algorithm in Appendix B, and may be computed only after the optimal G is obtained from (30). 
The following theorem states the necessary conditions of the optimization problem (30). 

Theorem 2 : 

Necessary conditions for G to be the solution of the optimal controller design problem (30) are: 



[A + BGM]Xj [A + BGM]* + DW p D* + (BGIi )(V p + E r )(BGIi )’ + BIj ^(Bl! )* - X! = 0 (31a) 


[A + BGM]*Ke[A + BGM] + [C + IoGM]*Q[C + IoGM] - Ke = 0 (31b) 

[A + BGM]*K 2 [A + BGM] + [C + loGM]*Q[C + IoGM] - K 2 + V x = 0 (31c) 

[A + BGM]K 3 [A + BGM]* -K 3 + V k =0 (3 Id) 

(IoQIo + B*K 2 B)G(MX]M* + Ii(V p + E Z )I*) + (IoQIo + B KeB)GMK 3 M + 

+ B*(K 2 AX I +K e AK 3 )M*=0 (31e) 


where V x has 4 subblocks as 


V x (i, j) = 0 i *2 or j *2 


V x (2, 2)=-i-trI k ^ 
sn c 


"c Ke(2, 2)[E- 1 ]-rw[E]- C0 ] 


i=l 




and. V k also has 4 subblocks as 


V k (i, j) = 0 i * 2 or j * 2 


V k (2, 2) 


_ _q_ 

sn c 


tr Xir 


"c [E 

I- 

i=l 


” 1 ]Lw[E]* co iX 1 (2, 2) 


where E is the matrix of eigenvectors of the matrix K e (2,2) (2,2). 


□ 


The proof of theorem 2 is given in Appendix A. 

Remark 1: The only terms in (31) which are affected by q are the two terms in (31c) and (3 Id) 
denoted by V x , V k . Hence setting P = °° gives q = 0, V k = 0, V x = 0, K 3 = 0, K 2 = K^. Hence, 
eqs. (31) reduce to the standard LQG design by setting |3 = In this case, the 11 block of (31a) 
reduces to the Kalman filter Riccati equation, and the 22 block of (31c) reduces to the control 
Riccati equation. 



Remark 2: We shall denote the controller satisfying (31) as the LQGpw controller to indicate 
that the LQGpw controller requires an additional step; the computation of T c from Appendix B. 

Now, we have following LQGpw controller design algorithm: 


The LQGpw Algorithm 


Step 1: Solve G from equations (31 a)-(31e). This gives the LQGpw controller. 


Step 2: Compute T c = Ux£^UkU t Il,V* by solving U x ,E x ,U k ,U t ,n„V( from (28b)-(28e), 

using the G obtained in Step 1. 

Step 3: G = 

□ 


I 

0 


'I 0 ' 

0 

T -1 

G 

0 T 

— c 


— c 




is the optimal LQGpw controller for implementation. 


Remark: A natural algorithm to suggest in Step 1 is as follows. Suppose one desires to design a 
LQGpw controller for 10 bit arithmetic. 

(i) Solve (31a)-(31e) for ft = «>, (hence, the standard LQG controller). 

(ii) On the next iteration set ft = 32 (or whatever gives a reasonably small number for 
V Vv 

(iii) Iterate by indexing ft. Change ft by no more than one bit on each iteration. This gives 
an "answer" in 32-10 = 22 iterations (but this manner of choosing step sizes in not 
guaranteed to be sufficient to yield the optimal answer). 

This is a "natural” homotopy method, since (3 is a natural choice for a homotopy parameter. 
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V. Computation Examples 


We consider an Euler Bernoulli beam modeled by its first 5 bending modes with 2 inputs 
and 2 outputs. The modal frequencies appear in TABLE 1. In discrete controller design, the 
discrete model is represented by the matrices { Ap,Bp,Cp,Dp,Mp,W p ,Vp } in equation (1). 
These matrices are given in Appendix C for a uniform sample time At = 0.018 sec. The LQG 
cost function is given by equation (7) with 

Q p = 0.991 R = 0.011 . 

The wordlength of the control computer is assumed to be 4 bits. Since the effects of D/A and 
A/D conversion errors on the control system simply modify the effects of system disturbance 
and measurement noise, we ignore these errors in the example. Both the standard LQG 
controller and the LQGpw controller are computed for the system. 



Frequency 

Damping Factor 

Mode 1 

3.4987e+00 

9.9994e-03 

Mode 2 

1.3995e+01 

2.1301e-02 

Mode 3 

3.1488e+01 

4.5600e-02 

Mode 4 

5.5979e+01 

8.0400e-02 

Mode 5 

8.7468e+01 

1.2530e-01 


TABLE 1. Frequencies and Damping Factors of the 
Euler-Bemoulli Beam Example 
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The standard LQG controller of course was designed without consideration of round-off 
oxtots ((3 = °°) and is labeled controller LQG in the TABLES. Controllers denoted LQGy. 
i = 1, • • • , 4 are the same as the LQG, but for a coordinate transformation on the controller after 
G is computed. The matrices ( A c ,B c ,C c ,D c } associated with the LQG T i controller are shown 
in Appendix C. In different coordinates Tj, TABLE 2 shows the finite wordlength contribution 
J e in the closed-loop system cost, using the standard LQG controller. In the optimal coordinate 
Tj (controller LQGti) the cost J e is about 500 times smaller than the cost in the original 
coordinate design (controller LQG). This improvement is equivalent to increasing the 

wordlength of the control computer by about 5 bits (5 = y log2500). The effect of 

computational errors J e in two commonly used coordinates, Normalized Observable Hessenberg 
Coordinates [Skelton 1988] and Phase Variable Coordinates, are also given in TABLE 2. The 
fact that Phase Variable Coordinates are bad for computation is consistent with other findings in 
filter synthesis [Williamson 1990]. The extreme high costs of the controller in a particular 
coordinate (LQGj4 ) in TABLE 2 serves only to demonstrate that the cost J e can become 
unbounded for some coordinates. The choice of coordinate T4 was rather arbitrary and will not 
be described or discussed further. 
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Controller 

Controller Coordinates 

Cost J e 

LQGn 

Optimal 

9.793 

LQG-n 

Normalized Obs. Hess. 

2.692 x 10 2 

LQG 

Plant Coordinates 

4.862 x 10 3 

LQG-n 

Phase Variable 

9.486 xlO 3 

LQGt4 

Coordinate "X" 

1.472 x 10 8 


TABLE 2. Standard LQG Controller in 
Different Coordinates 

The LQGpw controller was designed by the LQGpw algorithm given in Section 4. The 
controller matrices {A c ,B c ,C c ,D c } of this controller also appear in Appendix C. TABLE 3 
shows the computed costs of the standard LQG controller, the transformed LQG controller 
(LQGxi), and the LQGpw controller (The "LQGpw with coefficient error" will be discussed 
later). The costs for three different groups of excitations are computed in each case. The 
applicable disturbances for J, J y , and J u include plant disturbance w, sensor noise v, and finite 
wordlength error e. The applicable disturbance for J e , J ey , J eu is only e, and for J^, J^y, 
are only w p and v p (no finite wordlength effects). Hence, these sums apply to the various cost 
decompositions; J y is the output term of J (the total cost), J u is the control term in J, hence J = 
J y + J u . Jwvy is the output term of (the contribution of v p and w p in J), where J y = J^y + J ey 
and Jg Jey Jeu* J J wv Je* Jwvu is the control term of J ^ and Jy — Jwvu Jeu* As we can 


no 


















Disturbances 

Applied 

Costs 

LQG 

Controller 

LQGti 

Controller 

LQGpw 

Controller 

LQGpw 

with coeff. errors 

All v, w, 
and e 

J 

4.8827e+03 

3.0589e+01 

2.1207e+01 

2.4695e+01 

m 

2.8053e+03 

2.3458e+01 

2.0798e+01 

2.4232e+01 

Ju 

2.0774e+03 

7.1303e+00 

4.0941e-01 

4.631e-01 

e only 

Je 

4.8621e+03 

9.9302e+00 

2.0067e-01 

1.4071e-01 

Jey 

2.7850e+03 

3.1790e+00 

1.3841e-01 

1.0275e-01 

Jeu 

2.077 le+03 

6.7512e+00 

6.2267e-02 

3.7961e-02 

v and w 
only 

J wv 

2.0659e+01 

2.0659e+01 

2.1006e+01 

2.4554e+01 

Jwvy 

2.0279e+01 

2.0279e+01 

2.0659e+01 

2.0279e+01 

Jwvu 

3.7912e-01 

3.7912e-01 

3.4715e-01 

4.2514e-01 


TABLE 3. Evaluation of LQG Controllers in Plant Coordinates, Optimal 
Coordinate and of the LQGpw Controller 


see in the TABLE 3, even when the standard LQG controller is in its optimal coordinate 
(LQGti), the J e portion of the cost is still about 33% of the total cost (9.9302 compared to 
30.589). By using the new LQGpw controller design algorithm, we reduce the J e portion of the 
cost 50 times, compared to the LQGti controller and 24,110 times compared to the LQG 
controller. In the latter case, this is equivalent to increasing the wordlength of the control 
computer by about 7 bits, That is, controller LQGpw will give the same performance using 4 bit 
arithmetic that LQG gives using 1 1 bits. Furthermore this improvement in output performance 
is accompanied by a reduction in control effort RMS = V.4094 1 vs. RMS = \ 2077.4 . To 


m 


































conclude this point, we see that if both controllers use 4 bits, the difference in RMS output 
performance is an order of magnitude (V20.798 vs. V2805.3 ). This kind of improvement in 
performance can mean the difference between feasibility and infeasibility of some control 
missions. 

With the new controller, the round-off portion J e of the cost is only 0.85% of the total cost 
as opposed to 33% for LQG. Now let us discuss the cost Jwv, which would be the total cost if 
the closed-loop system was only excited by measurement noise v p and disturbance w p . That is, 
suppose the LQGpw controller was designed for 4 bits, but evaluated using infinite bits. These 
are the conditions of the standard LQG design, since there are no disturbances in the evaluation. 
Jwv of the LQGpw controller is a little higher than that of standard LQG controller. The output 
term of the cost is also a little higher and the control term a little lower. These indicate that the 
LQGpw controller is a little more conservative than the designed standard LQG controller. This 
compromise in nominal performance allows robustness to computational errors. Note in 
TABLE 3, that the quantities that are optimized by the theory (under the given conditions) are 
shaded. 

In the design of the LQGpw controller, the equations (31a) to (31e) were solved iteratively 
by a gradient method. The standard LQG controller in its optimal coordinate (LQG-pi ) was used 
as the initial controller design for starting the iterative process. Figs. 1-3 illustrate the 
convergence process for the LQGpw algorithm, plotting the total cost J, the wordlength cost J e , 
the the output J y and input J u performances, versus iteration. The optimal coordinate 
transformation played a crucial role in reducing the round-off errors (reducing the enror by 3-4 
orders of magnitude) as shown in Fig. 2. This was expected because the transformation was 
formulated in the optimization problem. The LQGpw controller was obtained after about 300 
iterative computations, but note from Figs. 1-3 that after 120 iterations one might have stopped 
with little loss. 
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Coefficient Errors 

In the introduction we promised some evaluation of the effects of coefficient errors. We 
argued that even though the LQGpw controller is optimized only for state quantization it 
performs well with coefficient quantization as well. To show this we introduced coefficient 
errors in the controller by using 4 bit precision instead of infinite precision in the controller 
coefficients. The key issue here is this. Quantization errors in the state degrades performance, 
but does not destabilize, since the effect of e is just a disturbance (note that all controllers in 
TABLES 1 and 2 are stable). Coefficient errors can easily destabilize. Figure 4 shows the closed 
loop pole locations using the standard LQG regulator (using infinite precision). The system is 
stable as marked by the x’s. When the controller coefficients are implemented using only 4 bit 
arithmatic, some poles as indicated by the o’s in Fig. 4, are outside the unit circle. Hence the 
standard LQG controller is unstable using a 4 bit control computer. 

Fig. 5 shows the improvement in the LQG controller by its optimal coordinate 
transformation before synthesis. This is the LQGji controller. The poles (o’s) are in improved 
locations compared to Fig. 4, but the closed loop system is still unstable. The coordinate 
transformation helped but not enough. Fig. 6 shows the LQGpw controller when controller 
coefficients are implemented using only 4 bits. The system is stable, confirming for this 
example improved robustness to controller coefficient errors, even though the controller has 
been optimized only for errors in controller state computation. The performance degradation in 
J, listed in the column "LQGpw with coefficient errors" in TABLE 3 is about 15% (compared to 
nominal performance in TABLE 3). 

Finally, we consider errors in both the plant and controller coefficients (due to 
quantization to 4 bits). These results are summarized in TABLE 4, where the modal damping in 
all modes is multiplied by parameter p. Hence p=l corresponds to the nominal plant in all of the 
prior discussion. The range for stability using the LQGpw controller is .729 <p< 1.23, 
demonstrating improved robustness over standard LQG controllers in the presence of errors in 
plant and controller coefficients. 



Damping Error Factor p 

LQG 

Controller 

LQGji 

Controller 

LQGpw 

Controller 

1.5242e+00 

unstable 

unstable 

unstable 

1.3717e+00 

unstable 

unstable 

unstable 

1.2346e+00 

unstable 

unstable 

STABLE 

1.1111 e-MX) 

unstable 

unstable 

STABLE 

1.0000e+00 

unstable (Fig 4) 

unstable (Fig 5) 

STABLE (Fig 6) 

9.0000e-01 

unstable 

unstable 

STABLE 

8.1000e-01 

unstable 

unstable 

STABLE 

7.2900e-01 

unstable 

unstable 

STABLE 

6.5610e-01 

unstable 

unstable 

unstable 

5.9049e-01 

unstable 

unstable 

. 

unstable 


TABLE 4. Robustness Controllers with respect to modal damping 
(4-Bit Wordlength Controllers) 











VI. Conclusion 


This paper solves the problem of designing an LQG controller to be optimal in the 
presence of finite wordlength effects (modeled as white noise sources whose variances are a 
function of computer wordlength). This new controller, denoted LQGpw, has two computational 
steps. First the gains are optimized, and then a special coordinate transformation must be 
applied to the controller. This transformation depends on the controller gains, so the 
transformation cannot be performed a priori. (Hence, there is no separation theorem.) The new 
LQGpw controller design algorithm reduces to the standard LQG controller when an infinite 
wordlength is used for the controller synthesis, so this is a natural extension of the LQG theory. 
It was shown both theoretically and by example that the choice of controller coordinates 
significantly influences the effects of computational errors on the control system and that there 
exists an optimal set of coordinates in which to do these computations. Since we have not 
obtained a closed form solution for the LQGpw problem, design of the LQGpw controller by this 
algorithm requires significant computation. Hence, the improvement of the new controller is 
achieved at the expense of extra computational effort in design. 
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Appendix A 


1 . Proof of Lemma 3 

a) Using the singular value decomposition of T = U,Il t V*, then the constraint equation (15b) 
becomes 


(V t nr 2 V l ”) ii = s for alii 


from Lemma 1, above equation is equivalent to 


tr(n t 2 ) = sn. 


i p • 


Now, let us study the cost yof (15a). Using the inequality 


rtAA*)* [,I(AB >’ 


tr(BB*) 


(32) 


(33) 


we have a lower bound on y 

y- {u t n 2 u*p} = tr{(n t u t *A/p)(n t u t *^)*) 

^ [tr((Tl t U*'>/p')(U* [^jp - ]" 1 )*] 2 _ (tr{n t )) 2 
* tr^h/PrWh/P]- 1 )*} trip- 1 ) 


Now, to prove that y is unbounded from above, we prove that for any large scalar m > 0, we 
have y(T) > m for some T. Let us choose a T having the following fl t : 

fl t = diag (flj) such that 

ni=n 2 = ••• =^-2 = 4" 

vs 


and 


H6 


-1 


VmtrCP- 1 ) 
V2mstr(P _1 ) -1 


, n. =Vmtr(P- 1 ) 


where m is so chosen that 


1 

m > T- 

2str(P _ ) 


Then 


tr(f[ t ) — — s ( n p — 2) + 

i=l Ili 


2mtr(P~ 1 ) - 1 ( 1 

mtr(P _1 ) mtr(P -1 ) 


= sn r 


Hence the chosen T satisfies the constraint (33). Now, we have: 


(tr(n t }) 2 
Y “ tr{P -1 ) 


"p - , 

(Zno 2 

i=l 

trlP- 1 } 


> 


(ri„,) 2 

— = m 

trlP" 1 ) 


we then conclude the proof of part a). The proof of b) and c) follows next. The lower 
bound and the matrix T are found by using following inequality: 


(trR) 2 < tr(QRQ* )tr(Q"‘ RQ _1 ) (35) 

the equality holds above when Q*Q = X, 2 I. 

Let us assume T = U t n t V*, P = U p n p U p , where Il t and n p are diagonal, U t , V t , U p 
are unitary matrices. Assume for the R and Q matrices in (35), 


R = u*u p n*u;u t (36) 

Q*Q = u:u p nJ i u;u t n 2 u;u p nJ i u;u t , (37) 


then 


(Q’Q)- 1 =u:u p n- 1/ *u;u l nr 2 u:u p n p ,/i u;u t . 

Hence, we have: / / 7 



tr(QRQ’) = tr(RQ‘Q) = trlCUrUpn^UjU.XUrUpn^UjU.nfurUpn^UjUt)] 

= tr[UpnpUju,n?u?] = tr[Pir*] =y 

tr(Q-‘ RQ- 1 ) = tr[R(Q* Qr 1 ] = uKurUpn^uJuoojrUpn^uJ^n^urUpn^u;^)] 

= trnj;u t nr 2 u:Up] = tr[nr 2 ] 

From equation (33), and the above equation we have the following: 

tr(Q-*RQ- 1 ) = tr[nr 2 ] = sn p 

Now, tr(R) — tr(U t UpFIp UpU t ) = tr(II^) = tr(Upn^Up) = tr'N^'. Substitute the above 
equalities back into inequality (35). We then have: [tr(VP)] 2 < sn p y, that is 

y > [tr(VP)] 2 

srip 

* 

Now, suppose the matrix T = U t n t V t yields the equality in (38) 
holds when Q*Q = X 2 I, then we have: 

u*u p n?u;u t n?u’u p n?u;u t = x 2 i , 

that is 

u«nfu* = x 2 u p n-*u; => u t nT 2 u t * = x 2 u p n*u; . ( 39 ) 

Hence 

-2_ Ppgg, 

X 2 

-2 

Substitute this n into equation (32) to obtain 

(v 1 o;u p F|-u;u 1 v;) U = s . 

= S “ P . hen “ = -J-«<nj?)« — n<VP). Now, substitute the above X 2 

Slip 

H8 


Then tr 


nVp 


X 2 


(38) 

. Since the equality in (35) 



into (39), to obtain 


_-2— * sn p u p n; /2 u; 


u,n t u t = 


tr(^) 


sn p >/p 

tr(VP) 


(40) 


Hence (38) yields the lower bound in (16a), and the matrix achieving this bound, shown by 
(40), must satisfy (17b). (17c) can be easily deduced from (15b). This concludes the proof. 

□ 


2. Proof of Lemma 4 


a) Proof of (29a): Since does not depend on K(p,q) for p i or q * i, we have: 


^k^rx(P’ Q) — 


3K(p, q) 


Jk.x =0 


Proof of (29b): We need following equality (e.g. Page 444 of Skelton [1988]) to prove the 
equation: 

MA] = rE _1 ]ith-rowA[E] ith _ col 


where \ is the ith eigenvalue of A, and E the eigenvector matrix of A. Now, we have by 
taking A = K(i, i)X(j, j) 

UK(i, i)X(j,j)] = [E-^/th-rowKG, i)X(j, jira/th-co, 


= tr[K(i, i)X(j, j)[E] ith ^ 


Hence from the differentiation rule 


StrAB _ 

dB 


we get 
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%vr "\ ~~ ^ ]/ih-row[E]/th-col^Cj’ j) 
oK(i, 1) 


Then, we have: 


^^k.x 

3K(i, i) 



^^go.oxgjM 

V^[K(i, i)X(j, j)] 

[E ]/ih-row[E]/th-col^O> j) 
VW i, i)X(j, j)] 


The proof of part b) follows in a similar manner 


□ 


3. Proof of Theorem 2 : Apply Lagrangian Multipliers K 2 , K 3 , then (30a)-(30c) leads to 
minimization of 


j = tr{Q([C + I 0 GM]Xj [C + I 0 GM]* + (IoG^XVp + E z )doGIi)‘)]} 

+ tr{K 2 ([A + BGM]Xj [A + BGM]* + DW p D* + (BGI^CVp +E I )(BGI 1 )* + (BIOE^BIO* 
- Xj)} + tr{K 3 ([A + BGM]*KJA + BGM] + [C + IoGM]*Q 

[C + I 0 GM]-KJ} + -2-(trI k ) 2 


sn c 


Then 
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dl 

5K 2 

dl 

3K 3 

JL 

dX x 

dl 

9Ke 


= [A + BGM]Xj [A + BGM]* + DW p D* + (BGIi)(V p + E z )(BGl!)* + BIjE^BIi )* 
= [A + BGM]*KJA + BGM] + [C + IoGM]*Q[C + IoGM] - K* = 0 
= [C + IoGM]*Q[C + I 0 GM] + [A + BGM]*K 2 [A + BGM] - K 2 + V Xi = 0 
= t A + BGM]K 3 [A + BGM]* - K 3 + V kj = 0 


-X! =0 


Applying Lemma 4 on the above two equations, we can obtain V Xi and V ki as stated in the 
theorem. This verifies (31a)-(31d). Now 


-1^- = 2l5QCX 1 M* +2 I;QIoGMX,M* +2i;QIoGI 2 (Vp + E Z )I* + 2B*K 2 AX!M‘ 
oG 

+ 2B*K 2 BGMX!M* +2B*K 2 BGIi(V p + E z )lI +2 B*KjAK 3 M* 

+ 2B*K 1 BGMK 3 M* + 2IoQCK 3 M* + 2i5qi<>gmk;m* =0 , 

but since IoQC = 0, then, 

= 2[IoQI 0 G(MX 1 M* + I a CV p + E z )lI) + B*(K 2 AX! + K 2 AK 3 )M* 
oG 

+ (B*K 1 B + IoQI 0 )GMK 3 M* +B*K 2 BG(MX!M* +I 1 (V p +E z )lI)] 

= 2[(lSQIo +B*K 2 B)G(MX)M* +I 1 (V p + E z )ll) + B*(K 2 AX 1 + K} AK 3 )M* 

+ (B*KiB + iSQI 0 )GMK 3 M*] = O . 

This verifies (3 le). 

□ 
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Appendix B 


We now present an algorithm (originally developed by Hwang [1977]) for solving (17b) 
and (17c) for one set of solutions of U t , TI t , V t (The solutions for U t , n t , V, are not unique). 
Let Vp" in (17b) be written in terms of its singular value decomposition 

^P=U p LpU; (41) 

where U p unitary, Lp diagonal. 

Algorithm (Solving U t , TT t , V, in (17b) and (17c)) 

I. Take: 


Ut = u p 


(42a) 



(42b) 

V sn p 



V, = V n _ 1 V n _ 2 ••• Vi • 

■ ■ V 2 Vj 

(42c) 


where V,, i = 1, • • • , n — 1 is computed as follows: 

II. Compute V j : Let 

I, 4nr 2 =diag(--- Ojj •••) (43a) 

Assume On and are two numbers such that one of them is bigger than s, the other is smaller 
than s. Then take V] as: 
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(3 column 


where 



ft 

0 ... 

0 

gl 

0 ... 

0 

Vi = 

0 

1 ... 

0 

0 

0 ... 

0 

0 


1 

0 



[3 row — > 

-gi 

0 ... 

0 

ft 

0 ... 

0 

0 

... 

... 


1 



0 





1 


fi = 


81 = 


o ip - 1 


Oip - On 
1-On 


°ip - o n 


% 


Vi 


(43b) 


(43c) 


(43d) 


Compute Vj : Let 




2-il ^i2 


(44a) 


where Z, e IR (l ^ satisfies the propeny [Z^ ^ = s, Z^ eIR (l I)x(n 1+1 ^ is a nonzero matrix, 
and Zj 3 can be written as 


£i3 = 


0 


0 

^nn 


Assume o^ and o ia are numbers such that one of them is bigger than s and the other is smaller 
than s. Then take Vj as 
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i column 


a column 


1 row 


Vi = 


a row 


1 0 0 

• . : 0 jO 

0 ... fi 0 ... 0 gi ... 0 

0 1 0 

o j jo 

0 1 0 

o ... gi 0 ... 0 fi ... 0 

0 j 0 j 

0 0 1 


Compute fj and gj as: 


fi = 


Pjq- 1 
CTia — Oii 




gi = 


f-gjj 
C^ia ~ ®ii 


1*4 


(44b) 


(44c) 


(44d) 

□ 


Computation of T c 

T c is formed as follows: = UjjL^UKUt^V* 

1) Compute the Covariance Matrix and Observability Grammian 

Ke = [A + BGM]*K e [A + BGM] + [C + I 0 GM]*Q[C + IoGM] 

X! = [A + BGM]Xj [A + BGM]* + DW p D* + (BGIjXVp + E z )(BGl!)* + BIjEuBIj 

Assume K e (2,2), Xj (2,2) to be (2,2) the subblocks of K e and X! (the controller 
subblocks). 
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2) Compute U x , L x , U k . 


These three matrices are computed by applying singular value decomposition on following 
matrices: 

X 1 (2,2) = u:i x U x 
Zi /2 U x K 2 (2,2)U;ri /2 =UkI k U k 

3) Compute U t> n t , V t . 

Let us replace P matrix in the algorithm of appendix B as 

P-diag [X i {K e (2,2)X 1 (2,2)}] 

Then we can compute U t> II,, V, by applying the algorithm on matrix P. 
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Appendix C 


DESIGN EXAMPLE OF ROUND-OFF LQG CONTROLLER 


Plant Model: 10th Order Euler-Bemoulli Beam 
Word-Length of the Assumed Computer: 4 bits 


1) The 10th Order Euler-Bemoulli Beam Model for Controller Design 


0 

0 

0 

0 

0 

0 

0 

0 

0.0095 

0.0937 


’ 0.9980 0.0179 0 0 0 0 0 0 0 

-0.2196 0.9968 0 0 0 0 0 0 0 

0 0 0.9687 0.0177 0 0 0 0 0 

0 0 -3.4620 0.9582 0 0 0 0 0 

0 0 0 0 0.8469 0.0166 0 0 0 

0 0 0 0 -16.4457 0.7993 0 0 0 

0 0 0 0 0 0 0.5594 0.0139 0 

0 0 0 0 0 0 -43.6477 0.4340 0 

000000 00 0.1138 

000000 00 -72.4045 
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0.0014 

0.0006 

0.1557 

0.0716 

-0.0004 

0.001 1 

-0.0480 

0.1257 

-0.0012 

0.0013 

-0.1299 

0.1440 

0.0007 

0.0012 

0.0720 

0.1164 

0.0007 

0.0007 

0.0588 

0.0588 


0.0014 

0.0006 

0.1557 

0.0716 

-0.0004 

0.001 1 

-0.0480 

0.1257 

-0.0012 

0.0013 

-0.1299 

0.1440 

0.0007 

0.0012 

0.0720 

0.1164 

0.0007 

0.0007 

0.0588 

0.0588 


0 7.8297 0 7.1091 0 -1.3744 0 -8.3569 
0 6.2128 0 -8.7875 0 6.2128 0 0 


0 -6.2128 

0 -6.2128 


0 7.8297 0 7.1091 0 -1.3744 0 -8.3569 0 

0 6.2128 0 -8.7875 0 6.2128 0 0 0 




1.0003e-03 

0 


0 

1.0003e-03 


2) Designed Regular LQG Controller in Optimal Coordinate LQG T i 
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A c = 


B c = 


C c = 


-0.4582 

-0.1633 

0.4144 

0.6040 

0.0849 

-0.5217 

0.4753 

-0.3503 

0.3326 

0.0383 

0.2946 

-0.1855 

1.5034 

-0.2726 

0.5293 

0.0908 

-0.0468 

-0.0574 

-0.4312 

0.1539 

’ 0.1894 

-0.2895 

-0.422 

0.0230 

-0.0296 

0.0941 

-0.0120 

-0.0024 

-0.0258 

0.0940 

-0.061 1 

0.0609 

-0.2200 

0.4919 

-0.0737 

0.2522 

0.0252 

-0.0076 

0.0737 

-0.0776 

-1.9370 

3.8601 

1.6850 

-3.2381 - 


-0.0133 

0.4587 

0.5622 

0.2226 

-0.5299 

-0.0.850 

-0.0095 

-0.0359 

-0.0709 

-0.0256 


-0.1836 

-0.4122 

-0.3257 

0.5105 

-0.1864 

-0.3095 

-0.2270 

-0.0617 

-0.0716 

0.0559 


0.1574 

- 0.0201 

0.3373 

-0.3084 

0.4324 

-0.2941 

-0.0416 

-0.3343 

-0.0416 

-0.1463 


-0.4386 

-0.0411 

0.2351 

0.0821 

0.3391 

-0.0605 

-0.4845 

-0.0787 

0.1318 

0.4745 


D c = 


-0.1054 

0.2748 

0.0665 

0.4446 

0.3306 

-0.7404 

-1.5704 

-1.0273 

0.5827 

-0.0777 


-0.2805 

0.1059 

0.1975 

0.1978 

0.2351 

0.0085 

-0.3867 

-0.1971 

-0.9215 

-0.3449 


0.2304 

-0.0786 

-0.1651 

-0.1382 

-0.1635 

0.1530 

-0.0236 

-0.0491 

-0.0746 

-0.9854 


-0.2815 

0.0379 

0.2658 

-0.0456 

-0.1155 

0.5389 

-0.4084 

0.4129 

0.2806 

-0.6735 


4.1659 3.4458 1.8923 -4.2436 -15.7358 -6.5380 3.7048 -5.3330 
•2.7357 -2.8624 -2.7744 5.2406 11.5365 3.9625 -2.8745 -0.1711 


3) LQGpw Controller from the LQGpw Algorithm of Section 4 
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i i i 


0.3501 0.4306 -0.2223 0.3078 -0.5350 0.1231 0.1595 -0.2003 -0.1024 0.1325' 

.2004 -0.2851 -0.2294 0.1810 0.6715 -0.4432 0.2756 0.1591 0.3525 -0.3974 

.2033 0.2556 -0.0197 -0.8326 -0.8293 0.0885 -0.0605 0.2870 -0.0571 0.1147 

.2973 -0.3621 0.6480 0.3770 -0.4095 0.4031 -0.2736 0.0125 0.0426 -0.0372 

0.0308 -0.2207 -0.5168 -0.3001 -0.9847 1.1705 -1.0703 0.7456 -0.0979 0.1920 

Ac= -0.1187 0.4836 0.0470 0.3655 0.2493 -1.0109 0.3516 0.5930 0.2744 -0.3872 

0.0089 -0.3363 0.0664 -0.0869 0.0085 -0.0712 -0.1936 0.113 0.3818 -0.5248 

-0.3341 0.2935 0.1055 0.1309 0.2251 -0.3631 -0.7912 -0.5655 -0.2610 0.3180 

0.0731 -0.0312 -0.0788 -0.1349 -0.4369 0.2594 -0.4096 -0.3895 0.7609 0.3237 

-0.1129 0.0070 0.0781 0.1679 0.4955 -0.3354 0.5865 0.4685 0.2460 0.6396 


-0.0134 0.0927 

0.0812 -0.1630 

0.0706 -0.3987 

0.2464 -0.6411 

0.4583 -1.0134 

Bc " -0.5942 1.0745 

0.2455 -0.2146 

0.1121 0.0815 

0.1013 -0.2475 

-0.1510 0.3465 


0.8861 -1.8997 3.8592 -0.3107 5.3072 -0.7395 0.6339 -1.6517 0.9202 -1.3734 

■1.4019 2.2532 -2.6576 0.1575 -3.3358 1.2007 -0.5179 0.2062 -0.9969 1.3884 



-0.4486e-04 -0.1328e-04 

-0.5913e-04 -0.1567e-04 
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